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The low temperature surface resistance R s of d-wave superconductors is calculated as function of 
frequency assuming normal state quasiparticle mean free paths £ in excess of the penetration depth. 
Results depend strongly on the geometric configuration. In the clean limit, two contributions to 
Ra with different temperature dependencies are identified: photon absorption by quasiparticles and 
pair breaking. The size of nonlocal corrections, which can be positive or negative depending on 
frequency decreases for given £ as the scattering phase shift <5n is increased. However, except in the 
unitarity limit <5n = 0.5 tt, nonlocal effects should be observable. 



I. INTRODUCTION 



With the onset of superconductivity, a new leagth 
scale, the BCS coherence length £ = vf/7tAo enters.Ellts 
size, rather than the size of the normal state quasiparti- 
cle mean free path £, relative to the penetration depth A 
determines whether or not the electromagnetic response 
of the condensate is local. The situation is different for 
thermally exited quasiparticles which are responsible for 
losses at frequencies below the pair breaking threshold. 
These quasiparticles behave very much like quasiparti- 
cles in the normal state once the change in the density of 
states and the resulting change in thermal occupation of 
quasiparticle states are taken into account .0 Hence one 
must expect to see nonlocal effects in the surface resis- 
tance of superconductors when I > A. At finite frequen- 
cies the quasiparticles contribute to the screening cur- 
rents so that there are some nonlocal corrections to the 
penetration depth of microwave holds as well. 

High-T c materials are highly anisotropic. They can be 
viewed as stacks of weakly coupled conducting planes. In 
this paper we shall only consider the in-plane conductiv- 
ity because it is only this component of the conductivity 
tensor which can be expected to show nonlocal effects. 
The in-plane coherence length typically is between 10 A 
and 20 A while the in-plane penetration depth is two or- 
ders of magnitude larger, so that hiflh-T c superconduc- 
tors are well into the London limitEJ, £ <C A. While 
the quasiparticle in-plane mean free path I is compara- 
ble to £ near the transition temperature T c , so that the 
corresponding scattering rate T = vp/2£ is of the or- 
der of T c , it is now generally accepted that £ increases 
substantially below T c .BlI In high quality single crystals, 
the disorder induced mean free path £ cl can exceed -the, 
London penetration depth by quite a wide margin.tri2'EZl 
In contrast to the normal state where a single param- 
eter suffices to characterize the elastic scattering, both 
the concentration of scattering centers and the strength 
of the individual scatterer or, alternatively, the normal 
state elastic scattering rate T and the scattering phase 



shifts have to be specified in the case of superconductors 
with unconventionalO order parameters. We shall show 
below for the case of s-wave scattering that the scattering 
phase shift 5n plays a very important role. 

It has been pointed out by Chang and Scalapinoi that 
nonlocal effects are absent when the conducting planes 
are parallel to the surface exposed to the microwave 
radiation (c-axis orientation) because the spatial varia- 
tion of the fields parallel to the surface occurs on length 
scales much larger than the penetration depth. Thus, for 
strictly 2 D quasiparticle motion, the scalar product -q 
vanishes. 

Zuccaro et aL0 argued that, at least for YBa2Cu307_5 
(YBCO), coherent coupling between planes is not negli- 
gible so that there exists a finite component Vp of the 
Fermi velocity parallel to the surface normal. Assuming 
an isotropic gap, these authors calculated the nonlocal 
corrections to the surface impedance for parameters such 
that £ <C A. These calculations show the expected in- 
crease in the surface resistance in the anomalous skin 
effect regime w < vp/\ = (£ /A)7rA <C A resulting 
from the direct absorption of photons by quasiparticles 
in pSJi energy- and momentum-conserving process. Li et 
aZ.til also invoked finite dispersion perpendicular to the 
planes to argue that nonlocal effects prevent the obser- 
vation of thje nonlinear Meifiner effect predicted by Yip 
and Sauls £3 While this assumption of a 3 D tight-binding 
Fermi surface appears quite reasonable for YBCO, one 
nevertheless expects v F <C f| b . Nonlocal effects will, 
therefore, be less significant for c-axis oriented samples 
than for samples whose surfaces contain the c-axis, like d- 
axis oriented films or thin needles of single crystals with 
the c-axis being the longest dimensions. Setting aside the 
problem of preparing high-quality a-axis oriented films, 
microwave measurements of the surface impedance of 
such films are not well suited to a search for nonlocal et 
fects because the typical current and field distributionsEl 
involve two vastly different components of the conductiv- 
ity tensor. Transmission experiments in the THz regime- 
using linearly polarized light, would be more feasible.t2l 



1 



For this reason we investigate the surface impedance for a 
rather broad range of frequencies. At very high frequen- 
cies, though, inelastic scattering would reduce the quasi- 
particle mean free paths even at low temperatures to such 
an extent that the local limit applies. An investigation 
of the electromagnetic response of a needle shaped sin- 
gle crystal to a parallel microwave magnetic field, which 
would involve the in-plane conductivity only, appears to 
be even more promising in view of the thick single crys- 
tals that can now be grown in BaZrC>3 cruciblesO 

If the pairing state in high-T c superconductors has 
nodes, which is at present a widely held belief, then a 
momentum dependent coherence length could be intro- 
duced. This would exceed the penetration depth for mo- 
mentum states in the nodal region so that high-T c materi- 
als are London superconductors only in the sense that for 
the majority of fc-states £ (fe).a^C A. This point was first 
raised by Kosztin and Leggett,Ej who found that the zero 
frequency clean limit penetration depth varies quadrati- 
cally with temperature at very low temperatures, rather 
than linear which is considered to be a hallmark of d-wave 
superconductivity. This change in the temperature de- 
pendence of A due to nonlocality is however easily masked 
by mean free path and finite frequency effects B 

Schopohl and Dolgo v0 have argued that a linear tem- 
perature dependence of the penetration depth would vi- 
olate the third law of thermodynamics. .—The nonlocal 
effects discovered by Kosztin and Leggettlij could recon- 
cile the d-wave model with this general thermodynamic 
argument but only for the geometry in which the. c-axis is 
parallel to the sample surface. Hirschfeld et al&3 pointed 
out that the electromagnetic response kernel contains, in 
addition to the vp ■ q term, a term q 2 /2m which leads 
to nonlocal corrections at some very low temperatures in 
any geometry. For q « 1/A <C kp this term can usually 
be neglected and this is what we shall do here because 
we are interested in microwave losses at low but finite 
temperatures. For strictly 2 D systems, one might won- 
der whether fluctuation effects would not be much more 
important than such very small nonlocal corrections. 

Another consequence of c?-wave pairing is the fact that 
pair-breaking is possible at any frequency. However, for 
this process to contribute to the transverse conductiv- 
ity, the finite photon momentum needs to be taken into 
account unless the required momentum transfer is pro- 
vided by some other scattering event occuring simultane- 
ously. Since we are interested here in low temperatures, 
we shall only consider disorder induced elastic scatter- 
ing. If only s-wave scattering is taken into account, it is 
easy enough to write down a general expression for the 
complex conductivity a(q,u>). In order to isolate the two 
contributions to the surface resistance which are due to 
nonlocality, we shall also consider the clean limit. 



II. THEORY, GENERAL CASE 

The current-current correlation function from which 
the transverse conductivity is derived according to 

2 

6 C Tt ~) 

= + (\jy>3y])(<l,u)j , (1) 

can be expressed very simply in terms ofj-jaormal and 
anomalous single particle Green's functionslla 

<bW,])(9,- m ) = 2T£j^§x 

Gu n (p) Gu> n -u m (p-q) + Fu> n (p) Fw n -v m ip-q) 

when there are no vertex corrections. This is the case for 
isotropic disorder induced scattering, which we will focus 
on here. Momentum dependent inelastic interactions, ca- 
pable of causing the formation of unconventional super- 
conducting .pair states, do require consideration of vertex 
corrections.Ej At low temperatures and low frequencies, 
however, the contributions of these interactions to the 
quasiparticle lifetimes are negligible. All the complica- 
tions resulting from quasiparticles in an unconventional 
superconductor being scattered off point defects, which 
were first discussed in the context of Heavy Fermion 
superconductors,!] affect only the single particle selfen- 
ergies discussed below. 

The momentum integral in Eq. (|^) is evaluated under 
the assumption that the main contribution comes from 
quasiparticle states near the Fermi surface. We consider 
purely two-dimensional conduction and simplify the 2D 
Fermi surface to a circle. The Fermi velocity vp and the 
density of states per spin at the Fermi level A^(0) = m/2it 
are combined into a single parameter, the plasma wave- 
length A p , according to 
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(3) 



where s is the average distance between conducting 
planes. In the clean limit, A p is identical to the zero 
temperature London penetration depth Al- With these 
assumptions one obtains for the conductivity after ana- 
lytic continuation of iv m to the real axis 
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where +(— ) indicates a positive (negative) infinitesimal 
imaginary part. M is an abbreviation for 
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The order parameter, which is assumed to have d x 2_ y 2- 
symmetry with the respect to the crystallographic axes, 
is represented as 



A(y?)=A (T) cos2(^-^ ), 



(7) 



where ip is the angle between the crystallographic (100)- 
axis and the surface normal. In the local limit, the ori- 
entation of the order parameter relative to the sample 
surface has no effect on the conductivity, provided one 
neglecte-,the suppression of the order parameter by the 
surface! 2 -!! which occurs on a length scale £ <C A p . The 
(^-integral can then be reduced to the interval [0, n/2]. 
In the nonlocal case with ip arbitrary such a reduction 
produces four different terms. For the exceptional cases 
ip = [(lOO)-oriented surface] and ip — 7r/4 [(110)- 
oriented surface] , these terms are pairwise equal. Further 
simplification is possible in the clean limit (see next sec- 
tion), as well as in the strong and weak scattering limit, 
because in these limits the selfenergy \ either vanishes 
or becomes independent of O, so that M depends only 
on [vpq cos ip} 2 . 

The dependence of M on q is very simple because we 
neglected q in the argument of the order parameter and 
the other self-energies. Anticipating |q| w 1/A P this is 
justified. If one wanted to keep, in addition to vpq cos tp, 
the term g 2 /2m,E3 such an approximation would be in- 
consistent. For q w fcp the response of a c£-wave super- 
conductor changes its character completely, because Iha. 
coherence factors would revert from case II to case lBt3 

Because of the restriction to s-wave scattering there 
are no self-energy corrections to the d-wave order pa- 
rameter. The remaining selfenergy corrections are to be 
determined from 
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(g ) is the energy- integrated normal Green's functions, 
averaged over the Fermi circle 



(fl (v.n±)> 



dip 



2tt 



! _ ttN(0)v 2 n 

N "imp . , 2 



1 + (irN(0)v) 2 kN(0) 



sin 2 <5 N 



(10) 



(11) 



is the elastic scattering rate in the normal state and 
<5n = tan _1 (7riV(0)v) is the scattering phase shift. 

From er(<7 = qe z , to) the surface impedance is calcu- 
lated assuming specular reflection 



dq 



7T J q 2 - iuj{i <j(q,uj) 




(12) 



We have not evaluated the formula applicable for diffuse 
surface scattering for lack of computer time. 



III. THEORY, CLEAN LIMIT 

Since we are primarily interested in the surface resis- 
tance, we shall calculate only the real part o\ of the com- 
plex conductivity En^^) . The imaginary part 02 should 
simply be given bya'til l/ui/igX 2 and this expectation is 
borne out by the numerical calculations based on the 
general formalism presented in the previous section. In- 
serting spectral representations for the Green's functions, 
evaluating the sum over Matsubara frequencies uj n and 
performing the analytic continuation with respect to iv m 
we obtain 
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the energy integral is easily done, yielding 
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+ a^n 2 - A 2 (ip) 



(17) 



So far, the evaluation ofj-sri is equivalent to that 
given by Mattis and BardeenEil who, at this point, take 
the extreme anomalous limit, i.e. they neglect the in- 
dependence of g a . If this approximation were used in 
the case of a d-wave superconductor with nodal lines 
parallel to the surface, the normal state conductivity 
<Ti = ^ JTT? ^vpQ ' a PP r0 P r i a te for the anomalous skin 
effect regime, would be obtained. Since HTC materials 
are essentially London superconductors, taking this limit 
is not justifiable. Instead, we shall find the zeros of g a (Q) 
and perform the 17-integral. In this way we cannot re- 
produce the normal state result because, for A — > 0, go- 
becomes either independent of fl or the 5-functions give 
a vanishing contribution. 

However, the same approach is used successfully in the 
theory of the electronic Raman response of high temper- 
ature superconductors, which involves a density-density 
correlation function with some vertex j(ip) specific to 
the Raman response in place of the current vertex p y /m. 
The important difference is that the coherence factor is 
case IjEJ i.e. A 2 in the numerator of Eqs. (p|JT^) has a 
different sign. One can then take the local limit q — > 
and evaluate the frequency integral using g~i(— w/2) = 
with the result 

- Im <[n, n]> (0, oj) = tanh (JjL} x 

x / ^ 7 2 (y)[2A(y,T)] 2 i 1 m ;IJ . (18) 

; 2 -[2A(^,T)] 2 



-■to, the Raman response 



This pair breaking contribution 

does account quite successfullycSEI for a range of exper- 
imental observations. The corresponding result for o\ is, 
of course, zero. This difference in the local clean limit 
explains why the absorption of low frequency photons 
by charge carriers in high temperature superconductors 
depends much more strongly on quasiparticle meanr-frec 
paths than the inelastic scattering of such photonsEj 
Solving <7±i(f2) = in the interval [— w/2, oo] leads to 



^o(f) = -77 + 77«FgCOS^/l : - 

2 2 y lo z — [v-p q cos ip\ 2 
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For 



- 9+i[-i] (fio(^)) to vanish, ^1 - _ ^^£ y]2 must 
be greater [less] than v F q cos tp so that either S [g+i] 
S [g~i] contributes. £lo((p) must be real, which requin 
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For costys > this is always fulfilled, irrespective of 
the magnitude and ^-dependence of the order parame- 
ter. rio(</3) is positive for such tp so that this represents 
the quasiparticle contribution to o~\ which, because of 
the Fermi function in (|l3|), becomes small at low tem- 
peratures. In the opposite case tto{tp) is negative. The 
inequality ( p0| ) then leads to the usual condition oj > 2A 
for pair breaking, when the order parameter is isotropic. 
For a d-wave superconductor, A(ip) has to be specified 
at this point. We introduce the abbreviations 



2A (T) 



t = cos ip, y 

For a (lOO)-surface we can then write 
A(y>) - A (T) [2t 2 - 1] 



v F q 



(21) 
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For this orientation of the d-wave order parameter the 
inequality (J20j) is fulfilled for t £ [i m i n , t max \ with 



max < 0, 



min < 1, 
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Such an interval exists only if the square root is real. 
For y < 1 <=> oj > 2Aq this is always the case. For 
y > 1, the interval is empty if x > Ay. For x < Ay, the 

square root is real provided x < Ay 



1 - i/i - l/y . For 

oj <C 2Ao this limits the g-integration in Eq. ( [l2|) to the 
interval [0, ^/2ui\ . 

Combining Eqs. @, @, and @ we obtain 

. 2(2A (T)) 2 



8 U/2oj -v F qj / dt + Q (v F q - u?j 



dt 



Vi-t 2 [f(n (t))~f(uj + n (t)) 
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(25) 
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where the first integral represents the pair breaking con- 
tribution while the second represents the contribution 
from quasiparticles. 

For a (llO)-surface, A(y>) = A (T) sin 2^, so that 
the factor \2t 2 — l] 2 in ( p5| ) has to be replaced by 
4 (l — < 2 ) i 2 . The limits of integration in the pair break- 
ing contributions also need to be changed. Instead of the 
interval [i m in,imax] around (p = 7r/4 we now have two 
intervals t g [0, t' min ] and t e [t' max , 1] near ip = and 
tp = tt/2 which can contribute to o\ 



given by expressions very similar to (23) and ( |24|) 
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IV. RESULTS AND DISCUSSION 



Parameters required as input for the numerical calcu- 
lations are the plasma wavelength A p , the in-plane Fermi 
velocity up, the transition temperature T c and the or- 
der parameter amplitude Ao(T = 0). The values chosen 
are typical of YBCO:Q T c = 90.5K, v F = 1.4 • 10 7 cm/s, 
Ap = 140 nm, and 2A (T = 0)/T c = 7.4. Because of 
(||), vp and A p are not entirely independent. The values 
given are consistent with an effective mass m rs 3 mo- 
T c sets the temperature scale but is irrelevant otherwise. 
The rather large value of Ao(T = 0) is deduced from the 
sharp drop in R s and A(T) observed in the vicinity of T c . 
If this were attributable to a non-BCS temperature de- 
pendence of Aq(T), a smaller value of Ao(T = 0) could 
be inferred which would increase all theoretical predic- 
tions for the low temperature surface resistance. 

A p and vp are the parameters which control the im- 
portance of nonlocal effects. We have performed some 
calculations with A p = lOOnm. Since in the local limit 
R s = 0.5 w 2 ^ApC7i is proportional to A p one has to scale 
pineal w jth \ i appreciate the differences in the nonlo- 
cal results. Because the effect of modest changes in A p 
are rather obvious, we shall not display these results. 

Fig. [l] shows results for R s at temperatures T = 0.13 T c 
and T = 0.05 T c for a disorder induced scattering rate 
r = 0.1 meV, which would correspond to a normal state 
mean free path of 460 nm. The scattering phase shift has 
been chosen as <5n = (Born approximation). Allow- 
ing for nonlocality yields peaks in R s at around 100 GHz 
which greatly exceed the results in the local limit. The 
peak heights decrease rapidly with decreasing tempera- 
ture, which indicates that these contributions to R s are 
due to direct photon absorption by thermally excited 
quasiparticles. This process ceases to be effective when 
u> ^> vpq. With q rs 1/A p we estimate that this "anoma- 
lous skin effect regime" should end at frequencies around 
v > i>f/27tA p rs 160 GHz in agreement with our numer- 
ical calculations. It is remarkable that here we have a 
range of frequencies in which R s is predicted to drop quite 
rapidly. 
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FIG. 1. Surface resistance as function of frequency for two 
temperatures, T/T c = 0.13 (dashed lines) and T/T c = 0.05 
(dot-dashed and dot-dot-dashed lines), for a d x 2_ y 2 order pa- 
rameter with its lobes perpendicular and parallel to the sur- 
face. The dotted line is the contribution to R s from pair 
breaking in the clean limit. The solid line represents the sum 
of this dotted line and the dot-dot-dashed line. 
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FIG. 2. Blow-up of the low frequency section of Fig. 
Pair breaking effects are negligible at these frequencies. 

At higher frequencies, R s increases again but is essen- 
tially temperature independent. This, one would suspect, 
is the contribution to R s from pair breaking. To check 
this, we calculated the pair breaking contribution in the 
clean limit from (f25|). The result is shown as dotted line. 
Taking the sum of this contribution and the result of 
the local approximation at T = 0.04 T c gives the solid 
line, which reproduces the result of the full nonlocal cal- 
culation remarkably well. The increase in R s found in 
the local approximation must also be attributed to pair 
breaking, made possible by the d-wave character of the 
pair state and the momentum uncertainty due to static 
disorder. 

While Fig. [l] shows the expected increase in R s , a com- 
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pletely different picture emerges when the same results 
are replotted for the small range of microwave frequencies 
relevant for applications (Fig. ||). Here we see that, de- 
pending on temperature, there exists a frequency below 
which R s in the local limit exceeds the results from the 
nonlocal calculation. This is a consequence of the partic- 
ular orientation of the order parameter we assumed. At 
these low frequencies the term vpq cos ip in (^|) empha- 
sizes the contribution from quasiparticles moving parallel 
to the surface. Since this is the direction in which the or- 
der parameter has its maximum, the number of occupied 
fc-states, that are effective in the absorption process, is 
reduced. The change in sign of the nonlocal correction to 
Rs with frequency would help to explain why theoretical 
predictions based on the local limit are higher than the 
low temperature data taken by Bonn et als on untwinned 
YBCO single crystals at 4.13 GHz and lower than those 
taken at 34.8 GHzO 

When R s becomes independent of frequency in the lo- 
cal limit, <7i must vary as 1/u 2 . In the normal state this 
would be the case for frequencies such that u/2T ^> 1, for 
parameters used here v 3> 50 GHz. In the superconduct- 
ing state, the effective scattering rate, which for given 
microwave frequency invoLves some temperature depen- 
dent frequency average,E§H depends very sensitively on 
<$n- For <5 N = 7r/2 it is usually found to be larger than T, 
while for <5n = it is much less than T [Ref. |5|, Figs. 14 
and 19]. Since <5n = was chosen to obtain the data 
shown in Figs. |l| and |2| it is not surprising that R s be- 
comes frequency independent at frequencies much lower 
than the 50 GHz estimated above. Since the surface re- 
sistance in c-axis oriented samples, where we believe the 
local limit to apply, does increase between 1 GHz and 
87 GHz, we must conclude that the scattering phase shift 
is closer to ir/2 than 0. 

The fact that the full nonlocal results for R s are lower 
at low frequencies than results obtained by taking the 
local limit does not depend in any essential way on the 
scattering rate or the scattering phase shift. This can be 
seen from Fig. || where both types of results, obtained 
for a range of scattering rates T and <$n = 0.4 n, are com- 
pared in a double logarithmic plot. While i? s depends 
very sensitively on T in the local limit, it is very nearly 
independent of scattering when nonlocality is taken into 
account, as one would expect in the extreme anomalous 
limit. Perhaps surprisingly then, and contrary to the 
behavior in the normal state, it is 7?' ocal that moves to- 
wards i?n° nlocal to close the gap between the two, when 
r increases. Unlike the high frequency regime shown in 
Fig. |l|, R s at these microwave frequencies is certainly not 
the sum of two independent contributions. The dotted 
line is the local result for the lower temperature already 
shown in Fig. [l]. Comparison with the other local curves 
indicates that in the superconducting state a small scat- 
tering phase shift is indeed more or less equivalent to a 
much reduced normal state scattering rate together with 
a large scattering phase shift . 
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FIG. 3. Double logarithmic plot of the local and nonlocal 
surface resistance for frequencies 1 GHz < v < 100 GHz and 
for different scattering rates F. The curves marked with num- 
bers are the results of the local approximations. The numbers 
themselves are the values of F in [meV] used in the calcula- 
tions. The four corresponding results of the fully nonlocal 
theory are those which attain the lowest values at 1 GHz. 
The dotted curve is the local result for F = O.lmeV and 
<5 N = 0. 
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FIG. 4. Surface resistance as function of the scattering 
rate for two different scattering phase shifts for fixed fre- 
quency v — 20 GHz. The inset shows Rs for <5n = 0.4 -k at 
v = 220 GHz. Local results are those that vanish for r — > 0. 
The temperature is T/T c = 0.04. 

The dependence of R s on T and <5n at a fixed fre- 
quency is further elucidated in the main frame of Fig. |J. 
In the unitarity limit = 0.5 n, nonlocal effects van- 
ish for r > O.lmeV. In this case the value of i? s can 
be calculated fiom the "universal" conductivity coo = 
l//i Ap7rAo(0).E3 In the local limit, R B must vanish for 
r — > 0. Before it does so, R S (T) goes through a max- 
imum, whose height, width and position depend sensi- 
tively on <5n- Because of this behavior one can fit the 
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low temperature surface resistance of high quality YBCO 
samples by choosing a small value of T which is, how- 
ever, still compatible with the expected structural disor- 
der, and then adjusting the scattering phase shiftl] With 
nonlocality taken into account, R s has a finite limit for 
r — * 0, as in the normal state. At the lower end of the 
microwave regime /jn° nlocal ^> i?' ocal holds only for values 
of r which would appear to be unreasonably small. As 
already shown in Fig. |2| at low frequencies one is more 
likely to find i?"° nl ° cal < i^al 

The inset shows -R s (r) at 220 GHz. At this frequency 
quasiparticles from the order parameter nodes can con- 
tribute to the photon absorption and pair breaking yields 
a sizeable contribution to R s so that nonlocal corrections 
are positive for all T, diminishing as T increases. 



0.3 



T = 0.005 meV (1 00) T/T c = 0.04 
5n = 0.4ti , - .----'I---'"' 




100 200 300 

v [GHz] 

FIG. 5. Surface resistance as function of frequency. Com- 
parison of results including a small amount of scattering ob- 
tained according to section |n] with results obtained by taking 
the clean limit. Dot-dot-dashed line: local limit. Dot-dashed 
line: nonlocality included. Solid line: quasiparticle contri- 
bution in the clean limit Dotted line: pair breaking 
contribution in the clean limit Dashed line: R B {v) in 
the clean limit. 

Fig. [| shows R B {v) for a very small scattering rate, 
but for finite 8^. In this figure one can clearly distin- 
guish the quasiparticle contribution (solid line) and the 
pair breaking contribution (dotted line) calculated from 
eq. (|25|). The sum of these two contributions, shown as 
dashed line, agrees well with the fully nonlocal calcula- 
tions according to section |[ These results are similar to 
those shown in Fig. |l| and demonstrate that an increase 
in 5n can be compensated by a reduction in T. 

Finally, we turn to the other extreme geometry in 
which the order parameter nodes are perpendicular and 
parallel to the surface [(llO)-surface, Fig. ||. In this 
case, nonlocal effects always lead to an increase in R s , 
but this is extremely small even if a very small value of 
r and <5n = is assumed. The corrections, visible in 
our calculations only at low frequencies, represent the 
contribution from quasiparticle states close to the order 



parameter nodes. Momentum conserving pair breaking 
processes are negligible in this geometry. This dramatic 
dependence of the nonlocal effects on the geometric con- 
figuration is easily understandable in terms of the clean 



limit formulae in section III 
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FIG. 6. Surface resistance as function of frequency for a 
different orientation, of the order parameter relative to the 
surface in the Born approximation for two different temper- 
atures (cf. Fig. [j]). For T/T c = 0.04 we also show results 
obtained by taking <5n = 0.4 n. 



V. CONCLUSIONS 



It appears to be possible that nonlocal effects can be 
observed in the surface resistance of high temperature su- 
perconductors, but only if the order parameter has nodes. 
For an order parameter with dz2_ y 2-symmetry only (100) 
or (010) surfaces can be expected to yield results signifi- 
cantly different from the local limit. Observation of such 
differences would provide direct evidence for the order 
parameter symmetry. Two contributions to the nonlocal 
response can be identified: absorption by quasiparticles 
and pair breaking. The latter contribution dominates 
in the THz regime, but since it is nearly frequency and 
temperature independent, it would be hard to identify. 

As in the anomalous skin effect in the normal state, 
the quasiparticle contribution to R s goes through a max- 
imum as the frequency is increased when <5n = 0. The 
peak height increases with the number of thermally occu- 
pied quasiparticle states until the temperature is so high 
that inelastic scattering severely limits the quasiparticle 
free mean paths. In a search for this effect it is therefore, 
not advisable to perform the experiments at very low 
temperatures. For the material parameter chosen the in- 
crease in R s over the local limit is largest for frequencies 
near 100 GHz. At the lower end of the microwave regime, 
nonlocality actually reduces R s below the local limit as 
a consequence of the anisotropy of the energy gap. 



7 



The size of nonlocal corrections to i? s in the supercon- 
ducting state depends sensitively on the scattering phase 
shift and not only on the disorder induced normal state 
scattering rate. If disorder were correctly described by 
the strong scattering limit, for nonlocal effects to be im- 
portant a degree of perfection in the Cu02-planes would 
be required that appears to be unattainable. When one 
moves away from this limit, the picture changes rapidly. 
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